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1. Introduction 

The findings in this paper are concerned with interval estimation of a location 
parameter 9 based on an observable X having cdf G{x — 9) with 9 > and 
symmetric densities g{x — 9), in other words for cases where there exists a 
lower bound constraint for 9 with the parameter space's minimal value taken 
to be without loss of generality. More specifically, we establish new analytical 
properties of the frequentist coverage 

C{9)^Pe{I^*{X)3 9), {9>0), (1) 

of the 100 X (1 — a)% Bayesian HPD confidence interval I^* {X) associated with 
the flat prior it* (9) = l[o,oo)(^)- 

With a known but conservative lower bound for the minimal coverage 
infe>oC(0) (see [9] for the normal case; [8] for the general unimodal case), we 
succeed here in establishing for logconcave pdf's g the better (for a < 1/3) lower 
bound 1 — ^ for inf6/>o C{9). We further estabhsh various new properties of the 
frequentist coverage C{9) for logconcave pdf's g, such as (i) argming>oC{9) = 
2do; (ii) 1 - ^ < infe>o C{9) < 1 - ^ + (Hi) sup,>o C{9) < 1 - f ; and 
(iv) C{9) decreases on {di,2do), increases on (2do,oo), where do = G~^{j^) 
and di — G~^(l — ^). Observe that result (ii) confirms the first order optimality 
of the lower bound 1 — as opposed for instance to the lower bound = 
1 — 2a + o{a) for small a. The main application of our results is undoubtedly 
for estimating a lower bounded normal mean with known variance, but the 
results are nevertheless more general and unified for symmetric densities with 
logconcave pdf's. For instance, an interesting illustration will be given for a 
Laplace model. 

Here is a glimpse of previous work pertaining to the credible interval Itt* {X) 
and its frequentist coverage C{9). For estimating a lower bounded normal mean 
9 (say 9 > 0), results due to (a) Roe and Woodroofe [9] (known variance), 
and (b) Zhang and Woodroofe [10] (unknown variance cr^) establish the lower 
bound for the frequentist coverage of the 100 x (1 — a)% HPD Bayesian 
confidence interval with respect to the priors: 7t{9) ~ l[o,oo)(^) in (s-); ^-nd 
no{9,a) = ^l[o.oo) (^)l(o.oo) ("■) in (b). More recently, Marchand and Strawder- 
man [8] established, for a more general setting with underlying symmetry, the 
validity of the lower bound for the frequentist coverage of the lOOx (1 — a)% 
Bayesian HPD credible interval derived from the truncation onto the restricted 
parameter space of the Haar right invariant prior. We refer to their paper for 
details and similar developments for non-symmetric settings (also see [11]; [12]). 
The important starting point to keep in mind from the results of Marchand and 
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Strawderman [s] is the applicability of the lower bound for symmetric and 
unimodal location models. 

As mentioned above, the analysis which we have carried out and which is 
presented herein was motivated by the conservative nature of the previously 
established lower bound and other unknown and potentially useful aspects 
of the coverage C (6) . Finally, the results are cast here with the backdrop of a 
flurry of recent activity and debate, which is focused on the choice of methods for 
setting confidence bounds for restricted parameters as witnessed by the works 
referred to above as well at that of Mandelkern [7], Feldman and Cousins [5], 
Efron [4], among others. 

2. Main results 

Our main results apply to location models X ~ g(x — 9) with densities g which 
are unimodal, symmetric about 0, and logconcave (i.e., log g is a concave function 
on its support). For such location families, the key assumption of logconcavity 
can equivalently be described as corresponding to those families with increasing 
monotone likelihood ratio densities, that is increases in x for all 9i,6n 

such that 9i > Oq. Alternatively, the assumption of logconcavity is connected to 
an increasing hazard (or failure) rate (sec Lemma 1). Before moving along, here 
is a useful checklist of notations and definitions used throughout. 

Checklist 

• 1 — a: credibility or posterior coverage or nominal frequentist coverage 
(as (0,1)) 

• g: probability density function (pdf) oi X — 6 

• G: cumulative distribution function (cdf) of X — 9 

• G(= 1 — G): survivor function 

• G^^: inverse cdf 

• A(-): hazard rate function given by \(z) = 

• TT*: the flat prior density truncated onto the parameter space [0, oo) given 

by^*(0) = l[o.oo)(^) 

• /tt* (X) = [l{X),u{X)]: the HPD Bayesian credible set of credibility 1 — a 
associated with tt* 

• a = lim^^_oo u{x) 

• G{9): the frequentist coverage at 9 of the confidence interval /^r* {X) given 
by Gi9)=Peih^{X)3 9) 

. di = G-Hl - f) 

• u~^{9) = M{x\9 < uix)} 

• 1-^9) = sup{a;|6l > l{x)} 



(these last two inverses being defined as such in order to facilitate expres- 
sions below for frequentist coverage) 
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We next collect some useful properties of logconcave densities (see [1], [2] for 
surveys) . Lemma 2 is a critical result which may well be new and of independent 
interest. 

Lemma 1. Let g be a unimodal, symmetric about and logconcave density on 
M. Then, 

(a) 

g{y+e) nondecreasing in 9 on [0, oo), for y > 0; 

(b) g^y^g-j nondecreasing in y on R, for 9 > 0; 

(c) g('J(gffe) ^nondecreasing in 9 on [0,oo), for y{-) a nonnegative and non- 
decreasing function taking values on [0,oo); 

(d) both the cdf G and survivor function G are log-concave on [0,oo); 

(e) the hazard rate A(-) is nondecreasing on R. 

Proof. Part (a) is obvious and is simply a consequence of unimodality. Parts (b) 
and (e) are log-concavity, and part (c) follows from (a) and (b). A proof of (d) 
can be found in [1] or [2]. □ 

Lemma 2. Let g be a unimodal, symmetric about and logconcave density on 
K. Then, for all z > 0, 

(a) 
(b) 



^ > S- - 1 ■ (2) 

g{2z) - 2G(z) ' ^ ^ 



^(^)>^-l. (3) 



G{2z) ~ 2G{z) 

Proof. Part (a) implies (b) given the increasing hazard rate property of part (e) 
of Lemma 1 , as 

X(2z) > X(z) =^ 3^ > 4^ for all z > . 
' G{2z) - g{2z) 



For part (a), begin by observing that for all z > 0: 

1 < =^ — =^ 



2G{z) -G(z)-G(2z) 
so that (2) holds as soon as, for all z > 0, 



9{z) A 5(22: 
or equivalently, 

' ,9{.y) giv + z) 



5(2) 5(2z) 
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Finally, apply part (b) of Lemma 1 inside the above integral to infer that |||y — 

^gj2zf — ^ ^'^'^'^ that y <E {0,z), which yields (4) and the proof 

of (2). □ 

We now recall previously established properties of J^. {X) and of its frequentist 
coverage before pursuing with the main analysis. 

Lemma 3 ([8], theorem 1). For X ^ gix ^ 9), 9 > 0, g unimodal, symmetric 
about 0, and the HPD credible set we have 

(a) h,{X) = [l{X),u{X)], withl{x) ^ {x-G-i(i + i^G(a;))}l(rf„,o,)(x) and 
u{x) = {a--G-i(aG(x-))}l(_oo,dol(^) + {^+G-H5 + ^G(a;))}l(d„.oo)(^); 

(b) C{e)>^,for alie>0; 

(c) G(0) = j^; 

(d) lime^ooGf^) = l-a. 

As shown in the following lemma, the lower bound l{-) is nondecreasing on 
(— oo, oo), while the upper bound u{-) is nondecreasing on (do, oo). Furthermore, 
we show how the logconcavity of g forces u(-) to be nondecreasing on (— oo, do) 
as well. (It is easily verified that Z(-) and u{-) are continuous functions on M.) 

Lemma 4. Consider It^-'{X) as given in Lemma 3. Then, 

(a) is nondecreasing on (—00,00); 

(b) u{-) is nondecreasing on {do, 00); 

(c) u{-) is nondecreasing on (—00, dp) whenever g is logconcave. 

Proof, (a) Since l{x) =0 for .t < do, we only need to look at the behaviour of 
l{x) for X > do(> 0). We have, for x > do, 

d 1 - a g{x) 



dx^^""^ ^ 2 5(G-i(i + i-^G(x)))- 

Notice, since l{x) > 0, that we must have x > G^^(i + i^G(a;)) > 0, from 
which we infer that g{x) < g{G^^{^ + ^-y^G{x))) given that g is unimodal with 
a mode at 0. 

(b) Follows directly with G and G^^ being nondecreasing.^ 

(c) Under the prior tt*. the posterior density of 6'|a; is given by 1[o,oo)(^)}; 
and u{x) is, for a; < do, the corresponding quantile of order 1 — a ([S], proofs of 
Theorem 1 or Lemma 5). Now the logconcavity of g implies that the family of 
posterior densities of 6*13::, with parameter x, has increasing monotone likelihood 
ratio in 9. Finally, the result follows since the quantiles w(x); x < do; of such 
families are nondecreasing in x. 

The increasing properties of and u(-), for logconcave densities g, permit 
us to translate the events of coverage, i.e., {/^-(X) 9 9}, as {u~^{9) < X < 
l~^{9)}, and leads to the following expressions and further properties of the 
frequentist coverage C{9). □ 



^actually, we have the stronger result that u{x) — x increases in x here for x > do 
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Lemma 5. Under the conditions of Lemma 3 with g logconcave, we have: 
(a) 

c{e) = <^ G{xi{e)) - G{xom ifO e (a, 2do] (5) 
[ G{xi{e)) - G{x2{e)) z/0e (2do,oo); 

where a = limi-_»_oo u(a;) (see checklist), and with xq{-), xi{-), and X2{-) 
being functions defined by the equations: 

ii) G{xo{e)) =aG{xo{e)+e) (0G(a,2do]); 
(m) 2G(a;i(6i))-l = {1 - a)G{xi{e) + 9) (6* > 0); 
{iii) l~2G{x2{9)) = {I - a)G{x2{e) + 9) (0 G [2do, oo)). 

(b) Furthermore, xq and xi are increasing, while X2 is decreasing, with 

lim xo[9) = ~oo,a;o((ii) = -di, a;o(2do) = -^o; 

9^0+ 

lim xi(9) ~ dn, lim xi(9) — di; 
61^0+ ■ e->oc 

X2{'2do) = —do, and lim X2{9) — —di. 

8—*QC 

Proof, (a) First, observe that the nondecreasing properties of Z(-) and u{-) es- 
tablished in Lemma 4 imply that l~^{-) and u~^{-) are equally nondecreasing. 
Along with Lemma 3, it also follows that l{x) varies continuously from to +oo 
as X G and that u{x) varies continuously from a to +00 as x £ 5R. Set 



xoi9) = u~\0) 

xi{e)^i-\e) 
X2{9) = u-\e) 

We hence obtain 



9 G (a,2do]; 
9 e [0,00); 
9 e [2do,oo). 



C{9) = Pg{l{X) <9< u(X)) 

= Pe{u-\9)<X <l-\9)) 

= G{l-\9)-9)-G{u-\9)-9), 

which is (5) with u^^(9) = —00 for 9 E [0,a]. Furthermore, it must be that: 

• l-^{9) - G-^{^ + ^G{l-\9)) ^ 9, for 9 > 0,^ G{l-^{9) - 9) = 
1 + i^G(ri(^))^'(*»); 

• u-\9) - G-'^{aG{u-\9)) = 9, for 9 e (a,2do),=» G{u~\9) - 9) = 
aG(w-i(6')) (i) 

• u-^{9) + G-i(i + ^Giu-^i9))) = 9, for 9 > 2^, =^ Giu-\9) -9) = 
\ + i^G(w-i(6')) ^ 

(b) The right hand sides of (i), (ii), and (iii) increase in given that both ^"-'^(0) 
and u^^{9') increase in 9. This implies that, in terms of 9 (> 0), xo(0), x\{&) in- 
crease, while X2(9) decreases. Now, observe that uid^) = 2do u~^{2do) = do, 
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and a;o(2(io)(= a;2(2(io)) = dQ—2dQ = —do. Similarly, u(0) = di ^ xo{di) ~ — di, 
and \imx-^-oou{x) ~ a ^ limg_>£j+ X(){6) = —oo. Now, making use of (ii) and 
(iii) and the limiting properties \\ms-,ool^^{Q) = oo, limg^tx) = C!0, we 

infer that \Ymg^ooXi{9) = di, and \v[a.s^ooX2{0) = — di. Finally, the prop- 
erty limg_>Q+ xi (0) = do follows directly from (ii), or again from part (a) of 
Lemma 3. □ 

Remark 1. Above, only the properties relative to xq, as well as those of the 
coverage C{9) for 9 < 2do require the logconcavity of g. Analogously, some 
results below (e.g.. Corollary 1) do not require the additional assumption of 
logconcavity. 

Remark 2. From (5) and the properties of xq, xi, and X2 of Lemma 5, it follows 
that C(-) is a continuous function on [0, oo), with the exception of a discontinuity 
at = a, and when a > 0. In this case, wc have limg^^- xo{9) = — oo, and 
lim9^(j+ xq{9) = u~^(a) — a, which will lead to a drop of G{u~^{a) — a) in the 
coverage at 9 = a. An interesting example of such a discontinuity occurs for a 
Laplace model (see Example 3, Figure 2, and Remark 3). 

Remark 3. It is pertinent here to discuss the behaviour of u{x) as a; — > — oo. 
In particular, we wish to single out cases where a > 0, which will imply that 
u{x) > 9 for any 9 e [0, a] (i.e., Itt'{x) does not underestimate 9 for such 0's, 
and coverage will occur as soon as underestimation does not occur; see (5)). As 
an example, consider a Laplace model with g{z) = G{z) ~ i e~^; for z > 0; and 
which leads to u{x) = — ln(a) for all a; < 0, hence a = — ln(a) > 0. Part (a) of 
Lemma 3 provides a way to verify this directly. Alternatively, observe that the 
posterior survivor function of 9 is given by 

P{9 > y\x) = - e-y; for y>0,x<0. 

Thus, for X < 0, the posterior distribution docs not vary and yields a constant 
I^.{x) = [0,-ln(a)]. 

Analogously, logconcave densities g with exponential tails will lead to a 
similar non-zero limit at — oo. A family of such densities, which will lead to 
u{x) — > — ln(a) as a; — > — oo, is given by g{z) = P(|z|)e~'^'; with P(-) nonde- 
creasing and logconcave on (0, oo), ^'(0+) < i-'(O), and — *■ as z — > oo. 
This may be verified by showing that the conditions on P force the density g 
to be logconcave, and that the posterior survivor function P{9 > y\x) converges 
(for y > 0) as above to e~y when x — oo. A simple example is given by 
P(z) = \i\z\ + 1), that is giz) = i(|z| + l)e-l^l. 

On the other hand, if h' is unbounded where h = — ln((7), then u{x) — ^ as 
X — > — oo. To prove this, it is sufficient to show that P{9 > y\x) — > as .t — > — oo 
for all y > 0. But notice that (for y > 0) 

hm Pi9 >y\x)^ hm = lim = 0, 



given that lim2;^oo{ h(y + x) — h{x) } > lim2;_>oo yh' {x) = oo. 
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Finally, we emphasize that the assumption of logconeavity is indeed required 
for the upper bound u(-) to increase on (—00, do). As for interesting counterex- 
amples, we point out that u{x) — > 00 as a; — > —00 for non-logconcave densities 
g such that lim2;_oo ^^g(t)^ ~ 1' foi' V > 0- This is the case for instance of 
Student densities with v > 1 degrees of freedom, as remarked upon previously 
by [10]. 

By virtue of Lemma 5, we now obtain a first set of new results for the coverage 

c{e), e>o. 

Corollary 1. For X ^ g{x — 6), > 0, g unimodal, symmetric about 0, the 
frequentist coverage C{9) of the HPD credible set (X) satisfies the properties: 

(a) C{6) increases in 6 for 9 > 2do; 

(b) Clff) <l - ^ for alie> 0; 

(c) C[e)>j^fore^[0,a]. 

Proof, (a) Immediate given Lemma 5's representation G{xi{6)) — G{x2{9)) for 
the coverage, and given that G{-) and a;i(-) are increasing on M, while X2{-) is 
decreasing on [2(io,oo). 

(b) It follows from Lemma 5 that C[e) < G{xi{e)) < G{di) = 1 - f , for aU 
6* > 0. 

(c) It follows from Lemma 5 that (7(6*) = G{xi{e)) > G(a;i(0)) = G{do) = 
for all 61 e [0,al. □ 

Corollary 2. For X ~ g{x — 9), 9 > Q, g unimodal, symmetric about 0, and 
logconcave, the frequentist coverage C{9) of the HPD credible set It,-* {X) satisfies 
the properties: 

(a) C{9) <l-a for all9> di; 

(h) C{9) >l-a for all9< G-\\ - ^) - d„ = d2 (say). 

Proof (a) It suffices to show that C{9) < 1 - a for e [di,2do) since C{9) < 
1 — a for all 9 > 2do, given part (d) of Lemma 3 and part (a) of Corollary 1. 
Observe that Lemma 5's increasing property of xq implies di < 2do given that 
xo{di) = —di < —do ~ xo{2da). Furthermore, the properties of Lemma 5 tell 
us that: xi{9) < di and xq{9) > —di for 9 G [di, 2do]. Finally as a consequence, 
we obtain from (5) with 9 G [di, 2do) ^ 9 > di ^ u{0) > a: 

C{9) = G{xi{9)) - G{xo{9)) < G{di) - G(-di) = 1 - a, for all 9 e [di,2do). 

(b) Take 6* > a as the result is already established for 9 < a, given that 1— a < 
j^. Now define 6*2 such that G{xq{92)) — It is the case that G{xq{9)) < 

for 9 < 92, given that xq{-) is a nondecreasing function. Furthermore, for 
9 < 92, we have 

C{9) > G{xi{9)) ~ G{xo{9)) 

I + a 1 + a L + a 1 + a 
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There remains to show that 62 = (^2 . But this is the case by definitions of Xq 
and da, with G-^{A) = -G-^{1 - A) for aU A e (0, 1), since 



= G(xo(02)) = aG(xo(02) + e2)^G 



XQ{e2) +92^02 = -do - G-' —— = d2. 

1 + a ' 



□ 



Remark 4. Observe that the analysis within the proof of part (b) of the previous 
corohary tells us also that 

sup Ci9) > C{d2) > G(.Ti id2)) - . (6) 

Example 1. For sake of illustration of the results obtained up to now. take a 
normal model with X\9 ^ N{9, 1), and 1 — a ~ 0.90. Evaluations give us di = 
G-i(0.95) ss 1.645, 2do = G-i(^) « 2.68, and = G-^{l-^)-dQ « 1.03. 
Results from Corollaries 1 and 2 thus imply that the frequentist coverage C{9) 
necessarily exceeds the nominal coverage 0.90 for 9 < d2 ^ 1.03; i.e., for values 
of the mean 9 that are within a little more than one standard deviation of the 
lower bound of the parameter space. On the other hand, the exact coverage 
falls below the nominal coverage and increases for 9 > di ^ 1.645. Furthermore, 
we may infer that: (i) G(0) = 0.90, (ii) supe>o G{9) < 0.95, and from (6) 
(Hi) supe>o G(6') > G(xi(d2)) - tx ~ G(a;i(1.03)) - ^ ~ 0.939. These 
above features, along with the increasing property of G{9) for 9 G [2(io,oo) 
(Corollary 1, part a), are illustrated in Figure 1. 

The next few results, culminating in Corollary 3, show that the minimum 
coverage is attained at 6* = 2do, and is bounded between 1— ^ and 1— ^ + iTq " 
The first of these results gives further information about the behaviour of C{9) 
for 9 <di. 

Lemma 6. For X ~ 17(2; — 9), 9 > 0, g unimodal, symmetric about 0, and 
logconcave, the frequentist coverage C{9) of the HPD credible set It^'^X) satisfies 
the property: 

C{9) > G(2do) for all 9 < di. 
Proof Take 9 e [Q,di]. From (5), we infer that C{9) > G(2do) as long as 

G{x,{9)) + G(xo(2do)) > G(xi(2do)) + G{xo{9)). (7) 

□ 

Lemma 5's properties of xq and xi imply that G{xi(9)) + G(a;o(2do)) > 
G(do) + G(-do) = 1, and G(a;i(2do)) + G{xo{9)) < G{di) + G(-di) = 1; which 
implies (7) and our desired result. 



E. Marchand et al./Frequentist coverage of Bayesian credible intervals 1037 

Lemma 7. For X ^ g[x — 9), 9 > Q, g unimodal, symmetric about 0, and 
logconcave, the frequentist coverage C{9) of the HPD credible set Ij^-lX) satisfies 
the properties: 

C[9) > C{2do) > 1 - Y' for all 9 > 2do. 

(ii) 

C(2do) < 1 - ^ 



2 \ + a 

Proof of (i). Since C(-) increases on [2do, oo) (Corollary 1, part a), it suffices to 
show that 

C{2d,)>l--. (8) 
From (5) and the other properties of Lemma 5, we obtain 

C(2do) = G(x-i(2do))-G(a;2(2do)) 

^ " " {G{xi{2da) + 2da) + G{x2{2d^) + 2do)} 



2 



i-ar„.„, 1 



1 + a 



From this, we obtain that condition (8) is equivalent to G(3do) > \-a'i ' ^'^'^ 
we conclude the proof by establishing the stronger property (given that do > 0) 

G 2do > Y ■ 

With the property G{dQ) = rewritten as a = i^G(do) ' infer that 



1 - a2 G{do) - l~a G{do) ~ 1 - 2G(do) 

Finally, the proof of (i) is complete since (3) implies the above chain of impli- 
cations, as well as (8). □ 

Proof of (ii). From (5) and the other properties of Lemma 5, we obtain in a 
straightforward manner: 

C(2do) = G(a;i(2do))-G(xo(2do)) 
< G{di)-G{-do) 



2 l+a 2 l+a 

The last piece of the analysis consists in showing that the frequentist coverage 
C{9) decreases on {di, 2do). The proof of the next result relies in part on several 
lemmas which are stated and proven in the Appendix. □ 
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Theorem 1. For X ^ g{x—9), 6 > 0, g unimodal, symmetric about 0, and log- 
concave, the frequentist coverage C{9) of the HPD credible set I^^-^lX) decreases 
on {di, 2do). 

Proof. Take 6 E {di, 2da). Make use (5) and part (c) of Corollary 4 to obtain 

C\e) < x[i9)gix,ie))~x',i9)gixoie)) 

^ (1 - a)g{da)g{dQ + di) 

- 2gido) - {1 - a)g{do + di) l-J^''' 

From this, we obtain that the property C'{d) < for 6* G {di, 2do) is equivalent 
to the inequality 

(1 - g{do) g{do + di) + a{l - a)g{di) g{do + di) < 2ag{do) g{di). 

Using the fact that di > dp > and the unimodality of g, we infer that the 
above condition is implied by either the condition 

g{di) > 1 - a. 



g{do + di) 2 a 
or, given part (b) of Lemma 1, by the stronger condition 

5(rfo) ^ 1 - a 



g{2do) - 2a ' 

Finally, the result follows with this very last inequality being equivalent to (2) 

Gjdo) 
l-G(do) 



given that a ~ ^j^,"} ^ ■ □ 



Corollary 3. For X ~ g{x — 9), 9 > Q, g unimodal, symmetric about 0, and 
logconcave, the frequentist coverage C{9) attains its minimum at 9 = 2da, and 
is bounded below by 1 — ^ . 

Proof. The result is a direct consequence of Theorem 1, Lemma 7, and Lemma 6. 

□ 



3. Examples and final comments 

Wc conclude with some comments and illustrative examples. We also refer to 
[3] and [(i] for further examples and illustrations. 

Remark 5. The new lower bound 1 — ^ for the frequentist coverage is an 
improvement over the existing lower bound for a < 1/3, and a significant 
improvement for relatively smaller a. For instance, with a nominal coverage 
of 1 — a = 0.90, the lower bounds are 0.85 and 0.81 respectively. Also, as 
alluded to in the introduction and as a consequence of Lemma 7, the bound 
1 — ^ is, for a < 1/3, fairly sharp especially for relatively smaller a given 
that M0>oC{9) < C{2do) < 1 - ^ + For instance with 1- a = 0.90, 
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□ 



d, 2 



Figure 1 : Frequentist coverage of !^(X) ; 
Normal model and 1 - a. = 0.90 



we obtain 0.85 < infe>o C{d) < 0.8590. There is some numerical evidence to 
support the apphcabihty of the lower bound 1 — ^ for some models that are 
not logconcave location models. This is perhaps the case notably for estimating 
a lower bounded normal mean with unknown variance, (a model associated with 
a pivotal Student distribution which is not logconcave), based on the reported 
numerical evaluations in [10]. 

Example 2. (Normal model; Example 1 continued) As a foUowup to Example 1, 
with X\e ^ N{e, 1), 6* > 0, and =l-a = 0.90, the coverage C{e) is presented in 
Figure 1. Notice that the coverage decreases on {di,2do) (Theorem 1), and the 
minimum is attained at 2dQ « 2.68 and is bounded above by 0.85 (Corollary 3). 
In fact, the exact minimal coverage C(2do) is about equal to 0.859, almost the 
same as the upper bound given above in Remark 5, or Lemma 7. 



Example 3. (Laplace model) Consider X\d ^ Laplace(^?) with g{z) = G{z) = 
\e^^ for z > 0, and G^^{t) = - log(2(l -t)) for t > 1/2. In view of these closed 
forms, explicit forms for /jr* (X) and its coverage C{9) are available. With further 
details provided in [A] or [* '•] , we obtain for instance from (5) , with a = — ln(a) 
(see Remark 2), do = G'^j^) = In(^), di = G-'{1 - f) = - ln(a) (= a): 



C{6) 



1 - 



2e8-(l- 



2(a- 



(l~a)e~ 



-a) 



if < 61 < - ln(a) 

if -ln(a)<^^<21n(iif) 
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□ 2 a 4 2do 6 8 10 12 

Figure 2 : Fre que litis t coverage of !Ij^*(X) ; 
Laplace model and 1-ot = 0.95 



Furthermore, we can show here that C{9) is precisely increasing for 9 € (0, di) 
and a < 1/3 (we aheady know it is decreasing on {di,2do), and increasing on 
(2do,oo)). This imphes here as well that supg>Q C(0) = C{di). A graph of 
this coverage is presented in Figure 2 for 1 — a = 0.95. We recognize several 
of the established features of C{9), such as those illustrated in Examples 1 
and 2. Observe the drop at a = -ln(0.05) « 2.996 of G{xo{a)) = G(xo(di)) = 
G{—di) = § = 0.025 (we made use here of (5) and Lemma 5's properties of 
Xq). Notice also that the minimum is attained at 2do = 2ln{:j^^^) w 4.70 
and equal to C(2do) ~ 0.92727 (using directly the expression above), which is 
very close to the upper bound 0.92738 given by Lemma 7. Finally, in this case, 
supg>Q C{9) = C{— ln(0.05)) « 0.9744; which of course is less (but only a little 
less) than Corollary I's upper bound 0.975. 

Finally, the findings in this paper does provide a sharper description of the 
frequentist coverage properties of the HPD credible interval I^{X) with the 
improved lower bound on the minimal coverage mitigating in favour of desirable 
features (of course, added to the fact that the interval I^{X) has exact credibility 
for a uniform prior on [0, oo)). However, one can turn around the argument to 
point out the non-conformity of the frequentist and nominal coverage of 
in the worse case scenario 9 = 2do. For instance, if 1 — a = 1 /3, Lemma 7 implies 
that infg>o C{9) is at most 7/12 (at least 1/2), in other words a departure of at 
least 1/12 between nominal and frequentist coverage a.t 9 = 2do. 
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4. Appendix 

All the results in this Appendix are called upon in Theorem 1 and are established 
under the same assumptions as those of Theorem 1, that is: X ~ g{x — 9), 9 > 0, 
g unimodal, symmetric about 0, and logconcave. 

Lemma 8. For 9 E [di, 2do], 

(a) the function g^^Q(0)+e) increases in 9; 

(b) the function g^xi{e)+0) "in^creases in 9, and is bounded below by -^j^fj^- 

Proof, (a) Follows from the unimodality of g and the properties of xq since 
Xo{9) < < xo{9) + 6* for 6* e [di, 2do]. 

(b) The increasing property follows from Lemma 1 (part c) and Lemma 5's 
increasing property of xi. With the addition of part (a) of Lemma 1. the lower 
bound is valid, since for € [di, 2do], 

9{xii9)) ^ g{xi{di)) ^ g{do) 



g{xii9) + 9) - g{xi{di)+di) " gido + di) ' 

□ 

Lemma 9. For 9 £ [di,2do], we have 

fn) t' (9) - (l-a)9(rri(9)+9) 
fh ) t' (f)\ — ag(xQ(e)+e) 

Proof. The expressions follow by differentiation of Lemma 5's implicit equations 
for xq and xi. □ 

Corollary 4. For 9 e [di,2do], 

(a) both x'i{9) and x'q{9) decrease in 9; 

(b) -'i(^) < .gJ„]Zt^Xit!-,d.) ' -o(^) > T^. 

(c) x',{9) g{x,{9)) < ,;;;„f^(tL^itC+:^V ""'^ ^o(^)3(-oW) > j^9{d^). 

Proof, (a) Follows directly from Lemma 8 applied to -p^jgj and -pjgj as given 
in Lemma 9. 

(b) The upper bound for x'l follows from Lemma 9 and Lemma 8's lower 
bound for -^f^^^^- Exploiting the fact that x'q is decreasing on [di,2do], 
Lemma 9, and the equality Xo(2do) = —do, the lower bound for x'q is obtained 
as: 

' (9)> ' (2d ) = ag{xo{2do) + 2^0) ^ ag{do) ^ a 

°' gixo{2do))-ag{xo{2do) + 2do) g{-do) ~ agido) l~a 

(c) Since xi (9) G [do, di] (Lemma 5), we have with g being unimodal and do > 0: 
g{xi{9)) < g{do) for aU 9 G [di, 2do] (and actually all > as well). The result 
then follows directly from (b). 
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Similarly, again making use of Lemma 5 and part (b), we have xo{9) G 
[-di,-do] for 9 e [di,2do], and x'„{0)g{xo{0)) > j^gi-di) = 1^5(^1), for 
all 6* e [di,2do]. □ 
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